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Abstract 

The effective interactions which violate a lepton flavor accompanied with neutrinos (nLFV) are 

> ' 

■^j- ' considered. Such a new physics effect is expected to be measured in future neutrino oscillation 

<N : 

experiments with long baseline. They are induced by radiative correction in the framework of the 

<N . 

. minimal supersymmetric standard model with right-handed neutrinos. We numerically evaluate 

in ■ 

the size of the couplings for nLFV interactions in this framework. The slepton mixing is not 
only the origin of the lepton flavor violation in the charged lepton sector (cLFV) but also that 
of the nLFV. We find that the nLFV couplings are strongly correlated with the corresponding 



X 



cLFV process, and they are constrained at O(10 5 ) times smaller than the standard four-Fermi 
couplings. 
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I. INTRODUCTION 



Numerous observations on neutrinos from the sun 
the reactor [101, and the accelerator 



Q, 3, 3, 0, 3, El , the atmosphere Q, 3, [J , 



111 ] suggest that neutrinos are massive and hence there 
are mixings in the lepton sector. This fact means that the standard model (SM) has to 
be extended so that the neutrino masses and the lepton mixings are introduced into the 
model. Lots of models to explain those experimental results have been proposed. Among 
them, a model with the seesaw mechanism^] has a promising attribute, in which tiny 
neutrino masses are naturally induced. Neutrino experiments also have revealed that the 
mixings in the lepton sector are much larger than those in the quark sector. This fact 
may imply that the lepton flavor number is strongly violated in the physics beyond the 
SM. Therefore, we can expect that the nature might exhibit sizable lepton flavor violation 
(LFV) and hence we could observe the remains of physics at the high energy scale. In the 
minimal supersymmetric standard model with heavy right-handed neutrinos (MSSMRN), in 
which the seesaw mechanism is realized, the LFV with charged leptons (cLFV) is expected 



to become large 
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15(. In this class of models, the renormalization effect due to the 



neutrino Yukawa couplings induces a significant size of off-diagonal elements of the slepton 
mass matrix, (rri-^a ( a P-> a i — e >/ i 7 r ); which are the seeds of the cLFV. Here, the 
flavor indices a and (3 should be understood to also indicate the mass eigenstates of the 
charged lepton fields. Concretely, the superpotential with the neutrino Yukawa couplings 
{fv)i a and the Majorana masses for the right-handed neutrinos M lJ includes the following 
terms, 

W D {f^WH^La + - 2 U' ! X,X r (1) 

where iVj, L a , and H u denote the chiral supermultiplet for the right-handed neutrinos, that 
for the lepton doublets, and that for the up-type Higgs field, and e is the anti-symmetric 
tensor for SU (2) l fundamental representation. The indices i and j are for the generation 
of right-handed neutrinos, which do not necessarily indicate the mass eigenstates. The 
renormalization group equation for the off-diagonal elements of the slepton mass matrix is 
given as l| , 

= 16^ + tif" m L + 2 ^lmlU + m 2 H jlf v + , (« + P), (2) 
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where m~ is the soft supersymmetry (SUSY) breaking mass matrix for the right-handed 
sneutrino, and m# is that for the up-type Higgs doublet. The matrix A v denotes the tri- 
linear scalar couplings corresponding to the first term in Eq. (fT|). Note that if the neutrino 
Yukawa couplings do not exist, then there is no LFV effect. We can approximately solve 
Eq. © as 



p _(6 + 2a|K , fi Mg 



(3) 



with a cutoff scale Mq and a typical mass scale for right-handed neutrinos Mr. Here, the 
universal soft SUSY breaking at Mg is assumed, and tuq is the parameter for sfermion 
masses and ao is that for scalar tri-linear couplings. In terms of the mass insertion method, 



X 




FIG. 1: One of the diagrams which contribute to the cLFV process lp — > l a ^f. This effect is 
approximately understood by the insertion of the LFV mass term (Am 2 ~)£ (a j3). 



we can see that the off-diagonal elements of the slepton mass matrix are the origins of 
cLFV. This is diagrammatically shown in Fig. From this diagram, it is obvious that the 
element (Am?)/ is relevant to the process lp —>■ l a ^j. In this class of models, the off-diagonal 
elements can become large so that the typical values of p redicted branching ratios are within 



a sensitivity reach of near future experiments [16. |17|,0]. Therefore, the search for the cLFV 
process is one of the promising ways to inspect the new physics effect beyond the SM. 

In this article, however, we investigate an alternative approach to explore the LFV, the 
search for the processes of the LFV with neutrinos (nLFV) at a long baseline (LBL) neutrino 
oscillation experiment. In the forthcoming experiments, the oscillation parameters such as 
the mixing angles and the squared mass differences are expected to be determined with 



high precision 
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2l|. Therefore, the measurement of nLFV effects might be possible. 
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The feasibility studies on the nLFV interaction search at future LBL experiments without 
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supernova neutrinos 
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30, fell. The 



assuming a specific model are made by Refs. 
current experimental bounds on nLFV interactions are given in Ref. [32j. The sensitivity 
for nLFV effects to solar neutrinos j^j, atmospheric neutrinos [3, the LSND results and 



36| has also been considered. It is pointed out that the nLFV signal is 



enhanced by the interference effect between the amplitude including the nLFV interactions 
and that of the standard oscillation (SO). 

We here investigate the nLFV interactions in the MSSMRN 1 . The origins of nLFV pro- 
cesses are the same as those of cLFV processes, which are the off-diagonal elements of the 
slepton mass matrix. They can become significant in this framework. In addition, there 
is an enhancement mechanism due to the interference effect. It can be expected that the 
detectable magnitude of nLFV effects is induced. We make the numerical calculations of 
the size of the nLFV couplings and show the correlation between the nLFV and the cLFV. 

In Sec. [Hi we recapitulate the model independent approach in the detection of nLFV 
effects at LBL neutrino oscillation experiments. We also show the way to parameterize 
the nLFV interactions. In Sec. we calculate these nLFV couplings in the MSSMRN 
and numerically evaluate the size of them under the universal soft SUSY breaking scenario, 
so-called the constrained MSSMRN. Here, we concentrate our attention on the nLFV inter- 



actions which are relevant in the oscillation channel — > v T . Finally, in Sec. IIVI we will 
give a summary. In Appendix we describe the model in order to make our notation clear, 
and in Appendix [B] we give formulae of the nLFV interactions. 



II. NLFV INTERACTION IN NEUTRINO OSCILLATION 

In this section, we explain how to parameterize the new physics effect with the model 
independent way|2y]. First, we note that in neutrino oscillation experiments we do not 
observe neutrinos themselves but do their products, corresponding charged leptons. Propa- 
gating neutrinos appear only in intermediate states. Therefore, the existence of the nLFV 
interactions suggests that there are some amplitudes whose initial and final states are the 
same as those of the SO which means the neutrino oscillation with SM interactions. 

1 The studies for the other models are done in e.g., Ref. [j^ . 
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To make the argument clear, we show an example. When we assume — > v T oscillation 
measurement at a neutrino factory experiment, all we can know are the facts, the decay 
of muons at a muon storage ring and the appearance of tau leptons in a detector which is 
located at a length L away from the source of the neutrino beam just after the time L/c, 
where c is speed of light. We interpret these events as the evidence of the SO, — > v T . The 
amplitude for this process A$o can be expressed by the product of the amplitudes for the 
sub-processes; 

Aso(^~ + 1 -> T ~ + F ) = A s{l^ -> v iJ p e e~)A v {v ii ^ v T )A d {v T d -> r~u), (4) 

where I (F) denotes all the other particles than a muon (tau lepton) in the initial (final) 
state, which can be measured in principle. In this example, J is d (a down-type quark) in 
Ad, and F means z/ e 2 and e~ in A s and u (a up-type quark) in Ad. The subscripts, s,p, 
and d, indicate "at the source of the neutrino beam" , "at the propagation process" , and 
"at the detection process" , respectively. Suppose that there is an effective four-Fermi nLFV 
interaction, 

£ e ff = A(e7 P P L /i)(^ T 7pP L i/ e ) + h.c, (5) 
we have the same signal through the other process than Eq. 

AiLfvGu~ + I —> t~ + F) = A s (fi~ — >■ v T v e e~)A v [y T — ° sc, > u T )A d {u T d -> r~u). (6) 

The external particles in Eq. © are completely the same as those in Eq. (jlj). Therefore, we 
can not distinguish the contributions from these two amplitudes. In quantum mechanics, 
we first sum up these amplitudes and next square the summation in order to obtain the 
transition rate. Therefore, an interference term arises between these two amplitudes for this 
process 3 : 

P{pT + I^t~ +F) = \A so \ 2 + 2Re[A* so A DLFY ] + \A nLFy \ 2 . (7) 

The term of the SO, the first term of the right-hand side, gives the leading contribution, 
and it is described by using the muon decay width T (= \A s (fi~ — > ^^ e e^)| 2 ) and the 

2 Precisely, it is one of the mass eigenstate of the neutrin oJSa . 

3 It is necessary to treat the neutrino as a wave packet |39j in the discussion on the coherence between 
these two amplitude. Here, we adopt usual treatment for the neutrino propagation, so that the neutrino 
propagation is described by the plane wave. 
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cross-section for the charged current interaction o (= \A&{y T d — > r u)\ 2 ) as 



\A 



so 



T x P v ^ Vr x a, 



(8) 



where P v -+ Vt is the oscillation probability for v» — > v T in the SO, which is defined as 



The second term which represents the interference between the amplitude 



of the SO Eq. (J3J) and that including the nLFV interaction Eq. 

A 



2ReL4* A lLFV ] = T x 2Re 



2V2G F A;{U " 



* v T )A v {v T 



x a, 



(9) 



where Gf is the Fermi constant. Note that the nLFV effect contributes to the oscillation 
probability not quadratically but linearly. Thus the effect can be enhanced and hence even 
if the nLFV coupling is small, it can contribute the oscillation probability significantly^. 

We now turn to the parametrization of effective couplings for nLFV interactions. As 
we have already shown, the amplitude for the neutrino oscillation process can be divided 
into three pieces, A s , A p , and Ad- First, we consider the decay process of parent particles, 
A s . Since all final states must be the same, nLFV interactions with (V — A)(V — A) type 
are important for the neutrino factory experimentQ. 

We can introduce the interference 
effect by treating the initial state of a propagating neutrino as the superposition of all flavor 
eigenstates. For the case of Eq. (J3J), we can take the initial neutrino state \u) as 



(10) 



where e^ T = A/(2 V / 2G F ). It can be generalized to the case of an initial neutrino with any 
flavor by using the source state notation which is introduced in Ref. j^] as 4 



/ I e s e s \ 



s i s 
^fie [it 



H) = (u')f\v a ), u s = 

C S _S 1 

We can include the total nLFV effect into the oscillation probability as 

2 



(11) 



PyS 



(12) 



The matrix U s is not necessarily unitary. 
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with the propagation Hamiltonian for the SO Hso] 



(S 



( 



SO)/3 



2E V 



(Ui 



MNSJ/3 







\ 



Am 2 21 



V 



MNS 



V 



o 



(13) 



where E v is the neutrino energy, a is the usual matter effect which is given as 2^/2G F n e E u by 
using the electron number density n e , ?7mns is the mixing matrix for the lepton sector, and 
Am^ (Am^) is the mass squared difference for the solar (atmospheric) neutrino oscillation. 
Next, we consider the propagation process, A p . The nLFV interactions modify the Hamil- 



tonian for neutrino propagation as 



24] 



nLFVj/3 



SO)l3 



2E U 



\ 



"ee "-"e/i 



a efi a fifi a fir 



\ci eT cl^ t (i TT j 



(14) 



where a a p is the extra matter effect due to nLFV interactions, which is defined as a a p = 
Sp= e du^VZ^/fGFripE,;, where n p is the number density for the particle p. Assuming the 
matter which consists of the same number of electron, neutron and proton, we can reduce 
it to 



l a(3 



2V2e^G F n e E p 



with 



€ a/3 = e af3 + 6e af3 + 6e a/3 ■ 



(15) 



(16) 



Note that to consider the magnitude of the matter effect, the type of the interaction is 
irrelevant since matter particles are at rest and hence the dependence on their chirality is 
averaged outj^. 

Then, we make a comment on nLFV interactions which affect a detection process, Ad- 
We can adopt a quite similar treatment to that at the source of the neutrino beam. In this 
article, we consider the case in which the nLFV interactions do not depend on the flavor of 
target-quark, which is almost the case for the so-called constrained MSSMRN. Therefore, 
we have the neutrino state for the detection process in the following form, 



\4) = [U d )p a \v a ), U d 
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J 



(17) 



Finally, the transition probability including the whole nLFV effects is given as 



P{fi +1 ->t +F)~Tx P v s a ^ v d x a, 



where 



(18) 



(19) 



III. NLFV INTERACTIONS IN THE MSSMRN 

In this section, we evaluate the effective couplings for nLFV interactions in the MSSMRN. 
The nLFV interactions are induced by the off-diagonal elements of the slepton mass matrix, 
which are similar to the cLFV interactions. We can naively estimate the size of the nLFV 
parameter e s a/3 from the diagram which is shown in Fig. l^asfilT] 



-a/3 



167T 



mi 



m s x \/ Br (^ Li)- 



(20) 



Here gi is the gauge coupling for SU(2)l and ms is the typical SUSY breaking scale. This 




FIG. 2: One of the diagrams which contribute to e^ Q . 



relation means that there is a correlation between two processes, VZ — > v a v e e~ and lp — > l a j. 
We here concentrate our attention on the nLFV associated with the tau lepton because 
that in \i-e sector is strongly constrained by corresponding cLFV processes. The current 
experimental bound on the branching ratio of r — ► /ry is 3.1 x 10~ 7 at 90% confidence 



level 



421 ]. and then this experimental bound constrains the nLFV coupling parameter e 



Mi- 



According to this naive estimation, the value of e* T may become as large as O(10 4 ). Such 
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size of nLFV interactions would be detected at a future LBL experiment such as the neutrino 

n 

factory [2 6|. 



A. Analytic Calculation of e s c ^ l ' d 

In this subsection, we explain the calculation to obtain e^,'™' d in detail and compare them 
with that of cLFV processes. The thorough results are given in Appendix [BJ 

For example, the effective Lagrangian relevant with the nLFV interactions which give 
potency signifies contribution to the oscillation „, in a neutrino factory is given 

as|26| 

+ 2V2G F e™?{v T Y?LV»)(Pl P p) (21) 

p=e,d,u 

+ 2V2G F ef lT (fYP L v,) (u lp P L d) . 

These effective couplings arise from penguin-type diagrams and from box-type diagrams as 
shown in Appendix [BJ The calculation of box diagrams is straightforward. It is almost 
the same as that of cLFV processes, e.g., fi-e conversion [3] except for the fact that only 
the (V — A)(V — A) type interactions are taken into account in the calculation of e s a/3 and 
elg. However, it is necessary to make an attentive calculation for penguin- type diagrams. 
In general, while the neutral and electromagnetic currents corresponding to Fig. ^take the 
following form, 

Pip - q) [i (2l(g 2 ) + <B(g 2 ) 75 } + {£(q 2 ) + V(q 2 )l 5 } Y + {(£(g 2 ) + d(q 2 )l 5 } <f] fp(p) 

+ h.c, (22) 

the charged current corresponding to Fig. El is decomposed to 

\p - q)P R {m'(q 2 )q u a^ + £'(g V + V} h(p) + h-C (23) 



pa/ 



Here, / denotes the charged lepton field I or the neutrino field u, p is the momentum for the 
incoming particle, and q is that for the gauge boson. All the coefficients, 21 (21'), 23, £ (€'), 
D, l£ (<£'), and are the functions of q 2 . In the limit q 2 — > 0, € and D for the electromagnetic 
current must vanish due to the gauge symmetry U(l) em . On the other hand, those for the 
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neutral and the charged current do not vanish because the corresponding gauge symmetry 
SU(2)l is broken. In other words, the Lorentz structure of the nLFV interaction which 
is shown in Fig. |2] is different from that of the cLFV process lp — > with a real photon 
emission. The former one is dominated by the vector exchange interaction (£'(0)), 5 while the 
latter is done by a di-pole type interaction (21(0) + 53(0)75). Therefore, they are definitely 
correlated each other but are not the same functions. 




FIG. 3: One example of the diagrams in which we have to take into account the off-shellness for 
the neutrino propagation. In the neutrino propagator which is pointed out by the oval, we must 
treat the neutrino as if it were a massless particle. 



In the calculation of nLFV diagrams, we have to pay attention to the following two facts: 
(i) we regard neutrinos as highly off-shell particle, and (ii) we must avoid counting one 
contribution twice. First, we explain the reason why neutrinos behave as highly off-shell 
particles, and hence it is necessary to treat that neutrinos are massless in diagrams for 
effective nLFV interactions. In Fig. ® fro,n the viewpoint of the uncertainty principle, the 
condition for neutrino oscillation to occur is described asj 

6(p-q)~— > , 24 

ox L p — q 

where Sx is the uncertainty of the position and the time in which neutrinos are produced, 
5(p — q) is the uncertainty of the energy- momentum of the outgoing neutrino, and Am 2 is 
the neutrino mass squared difference corresponding to the LBL oscillation experiment. The 
inequality results from the fact that the neutrino production position must be determined 
much more accurately than the baseline length L. The fact that the neutrino oscillation 



5 Since the effective four-Fermi couplings are induced by the exchange of the massive gauge bosons Z and 
W and q 2 is of 0(m/), essentially we can put q 2 = 0. 
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phenomena are observed in the LBL experiment requires that the equality at the most 
right-hand side should be satisfied. Thus, the uncertainty of the squared momentum must 
conform to the following relation 



Hip ~ <?) 2 } > Am 2 



then ((p — q) 2 ) > m 2 



(25) 



where ({p — q) 2 ) denotes the average of (p — q) 2 . This inequation shows that the average of 
the neutrino momentum is much larger than its mass, and it follows from this that neutrinos 
are generally highly off-shell fields in oscillation experiments. This also means that all the 
diagrams for the nLFV interactions include the off-shell neutrino as external lines. In the 
exact meaning of the field theory, we do not have the method to calculate diagrams with off- 
shell external legs. However, we can evaluate such diagrams by treating as if the neutrinos 
were massless. For the practical purpose, we make the following replacement which we refer 
as the off-shell prescription in calculations of the type of diagram shown in Fig. El 6 



Next, we turn to the problem of the double counting and show the way to solve it. First, 
we note that we must calculate the process pT + / — > t~ + F with the method of the field 
theory, so that we must not calculate the nLFV effect for each stage because we cannot 
observe the each stage separately. If we calculate the diagrams shown in Figs. IU-(a) for 
e* T andEJ-(b) for e™ T , then it means that we count twice the diagram of Fig. IU-(c) into the 
calculation of the process p~ + / — ► r _ + F. In order to avoid doubly counting, we have to 
get rid of one of them. For example, we should not include the contribution of the diagram of 
Fig. IU-(b) in . A similar situation occurs among e^ T and e™. The penguin contribution to 
e^ T is essentially given by the complex conjugate of that to the corresponding e% . However, 
we must eliminate the contribution from diagrams similar to Fig. EJ-(a) from e^ T . 

Finally, we should notice another thing that is important for the double counting problem; 
Which stage (source or matter) does the contribution of Fig. IU-(c) belong to? For example, 
instead of removing the contribution of Fig. IU-(b) in e™ , we can eliminate the contribution 

6 The mass parameter m Va for the flavor eigenstate[3^ appears in Eq. I|26|l . However, we finally neglect it 
in our off-shell prescription. 




(off-shell prescription) , 



(for usual on-shell particle case). 



(26) 
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(c) The whole of the diagram for the process fi + I — > r + F. 

FIG. 4: Schematic explanation for the double counting problem. If we count the diagram-(a) into 
e^ T and we also do the diagram-(b) into e™ T , then we count the diagram-(c) twice. 

of Fig. IU (a) in e* T . For this ambiguity, we adopt the way to divide the diagram into each 
stage so that epsilon parameters in each stage will disappear in the limit where mg — > oo 
after using the off-shell prescription. The SU(2)l symmetry is recovered and then e s j™ ,d 
should disappear in the large SUSY scale limit. This is shown analytically. Note that each 
diagram gives rather large contribution and stays almost constant in the large SUSY scale 
limit. The cancellation among the diagrams is highly nontrivial. Therefore, by checking the 
cancellation among the diagrams, we can be confident about the legitimacy of our treatment 
for the internal neutrino lines. 

B. Numerical Study 

A numerical study to evaluate the epsilon parameters is necessary in order to make it clear 
whether our naive estimation Eq. (|2(Jj) is correct. We here assume the universal soft SUSY 
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breaking at Mq and adopt the scenario of the radiative electroweak symmetry breaking |44j. 
The detail of calculations is shown in Appendices. We scan the values for the soft SUSY 
breaking parameters mo and M\/2 at the range of 100 - 1,000 GeV, and also scan the elements 
of the neutrino Yukawa matrix. We here take the normal hierarchy for the neutrino mass 
matrix. The scatter plots for size of e s , e™ T and e d are presented in Fig. |5J All of them 




Br(r — > fij) Br(r — ► fij) Br(r — ► /ry) 



FIG. 5: Scatter plots for size of nLFV parameters e^ r , e™ T and e^ r . Dots are for tan/3 = 10 and 
triangles are for tan (3 = 3. We here fix the soft SUSY breaking parameters as ao = and [i > 0. 

take similar behavior. As expected from Eq. (12011 . the nLFV parameters correlate with 
the branching ratio of the process r — > /ry. However, the size of the parameters is much 
smaller than that of the naive estimation. This fact can be understood by the cancellation 
among the diagrams which contribute to the nLFV interaction. In the SU(2)l symmetric 
limit, the diagrams for the penguin contribution to the nLFV interaction must cancel each 
other out because of the gauge symmetry. Since the diagrams for nLFV interactions are 
induced at the mg scale, in the limit where mg ^> mz, the SU(2)l symmetry is assumed 
to be recovered. Even in the case where mg ~ (9(100) GeV, the cancellation is rather 
significant. Thus, our naive estimation must be modified so that the additional suppression 
factor mz/ms is introduced. It may be worth pointing out that the epsilon parameters do 
not strongly depend on the value of tan j3 unlike the branching ratio of the cLFV process. 
This fact arises from the difference in the structure of the chirality in each process. The 
process r — > /ry is dominated by the diagrams including the left-right mixing of the slepton 
which is proportional to tan (3. The nLFV processes do not pick up such the left-right mixing 
term because a chirality flip is not necessary. Therefore, nLFV search may be advantageous 
in the case where tan (3 is relatively small. We also note that it is obvious from Eq. that 
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both nLFV and cLFV are enhanced when ao takes a large value. 
IV. SUMMARY 

We summarize our study and give some discussions. It is known that magnitude of nLFV 
couplings can become large enough to be detected at future LBL experiments within model- 
independent approach. We here considered the nLFV interactions in the MSSMRN under 
the universal soft SUSY breaking scenario which is one of the most promising candidates 
for the physics beyond the SM. 

We find that in this scenario the nLFV couplings cannot be significant, and hence it is 
quite difficult to observe these effects at future oscillation experiments. The reason why they 
are strongly suppressed is that the SU(2)l gauge symmetry is approximately maintained. All 
the particles and the interactions which can generate nLFV interactions respect the SU(2) l 
symmetry in the limit mj — > oo. Although each diagram to contribute to nLFV interactions 
can become large, a brilliant cancellation among those diagrams occurs. Therefore, the 
penguin contributions are strongly suppressed. 

We adopt the approximation for the calculation of penguin diagrams which is explained 
in Sec. IIII Al In order to confirm the validity of this approximation, we need to make the 
calculation for the process shown in Fig. IU-(c) by using the method of the field theory. 
However, the calculation which we adopt here must be reliable in the sense of the field 
theory because the consistencies, i.e., the recovery of SU{2) L gauge symmetry, are obviously 
maintained in our calculation. 

Finally, we mention our future work. Since the decay process of a muon differs from 
that of a pion, we might expect that sizable new physics effect exists only in the decay 
of pion. Thus, it is necessary to investigate nLFV effects in the MSSMRN at superbeam 
experiments individually. Furthermore, because within the model independent approach, 
epsilon parameters can be still significant and there are lots of the other models than the 
constrained MSSMRN which can explain the neutrino masses and the lepton mixings, we 
need to examine such possibilities. 
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APPENDIX A: MODEL 

We basically follow the notation which is adopted in Ref. jl^ . However, we pay attention 
to the fact that the mixing matrices to diagonalize the mass matrices for sferimons, chargino 
and neutralino are complex matrices, in general. 

The mass matrices for charged slepton, down-type squark and up-type squark are 

= (fl fn) (Mj): |7I= ltt In) I M b {M }/ 1171' ( A1 ) 




Mj Mj 

JLR JRR 




where 



ILL 



m \ + (#>) 75 Sin2 P + m ^ COs2 ^ { (a) - ^ sin 2 9 W ) , (/ = «), 
m l + (#>) cos2 /5 + ^1 cos 2/3 { (-1) - QL sin 2 6 W ) , (/ = J, d), 



(A2) 



M2 = 4 +(///}) ^ sin 2 /5- m| cos 2/3 {-QL sin 2 ^}, (/ = «), 
/M ™^ + (///}) ^cos 2 /3-m|cos2/3{-g/ m sm 2 ^} , (/ = l,d), 

M _ = i-A f ^sm(3-f,f f ^ 2 shx(3cot(3, (/ = «), ^ 
I A/^cos/^-yu/z-^cos/^tan/?, (f = l,d). 

Here, the indices a and /3 are for interaction eigenstates for thier superpartner fermion fields. 
We take the basis where the mass matrix for the charged lepton field is diagonal, so that 
the index for chraged lepton indicates its interaction eigenstate and its mass eigenstate, 
simultaneously. The unitary matrix Uf is defined as 

{Uf) x a {Mj)J{U}h Y = diag(m^ )S X Y . (A5) 

The relations between the mass eigenstates and the interaction eigenstates are 

fx = (U } ) x a f La + (U f ) x a+3 f Ra , (A6) 

f La = (Up a X f x , f Ra = (U}) a+3 X f x . (A7) 
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The sneutrino mass term is also given as 



where 



-Co = ** a (M?) a %, (A8) 



(Ml) =m 2 L + mlm u + m| cos 2(3 [ - ) , (A9) 



2 

where m u is the neutrino mass matrix which is induced by the seesaw mechanism 



m v = f u M- l f u V -smf3. (A10) 
The unitary matrix Up is defined as 

(^)x Q (M?)/(C/t)/ = d^g(ml x )5 x Y . (All) 
The relations between the mass eigenstates and the interaction eigenstates are 

VX = {Uv) X a VLa, "La = (^)/^. (A12) 

The chargino mass term in the 2-spinor representation is 



2 



M 2 \f2rriw cos (5 



h d J ^ " ' \ y/2mw sin j3 \x 



The diagonalization is done by unitary matrices Ur and Ul as 



)l 


(*-\ 







(A13) 



{U R ) A \M c )t{Vl)i B = diag(M- I )^ B , (A14) 
The relations between the interaction eigenstates and the mass eigenstates are 

{ZA)a = (U L U | ^" )8 | , (x +A ) a = (U R )\ ( { f ] a | , (A15) 

( h d)aJ i \{K) 

where a (a) denotes index for the Lorentz spinor for 2 (2) under SL(2,C). We here adopt 



the same rule as that in Ref. 



The 4-spinors for mass eigenstates are constructed as 



XI = I - )a ) > X +A =( { "l A A )a | , (A16) 



1(3 



and then those for interaction eigenstates are 
W~ = I 




H 




(w~ 



tm r i+ A} 



R) 

I) 1 A 



'x~ A ) 



(K)a \ = I {Ui) 2 \x- A ) a \ #+ = / (K)a 

(kr) W R ) 2 A (^A) & ) , \(h- d )\ 

The neutralino mass term in 2-spinor representation is 



T\2 (z,+A 



R) A' 



X 
X 



(A17) 



w° 



K 



w° K hi 



(A18) 



x 



/ Mi 





—mz sin 6*^ cos /3 


mz sin ^ sin (3 ^ 









M 2 


mz cos 6*vk cos /3 


— cos 9w sin /3 




w° 


—m z sin 6 W cos (3 


cos 9w cos /3 





-/X 




k 


\ m z sin 6w sin j3 


— m z cos 0jy sin /3 


-ft 


o ) 







The unitary matrix Un is defined as 



(^) A ,(^)^1)/ = diag(il%)<5 



AB 



(A19) 



The relations between the interaction eigenstates and the mass eigenstates are 



7,0 



X 



N A 



K 

\k) 



Using this 2-spinor x , the 4-Majorana spinor can be constructed as 



(A20) 



PlXa + PrX 



~0A 




(A21) 



The interaction eigenstates are 




TM t~0 A \a 



N) A 



( (t5°)J 


l-( 
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^0 = I ^d) a 



{U\,h A {x\) a 



(4)/(4)a 



The Lagrangian for gaugino-sfermion-fermion interactions is described as 



(A22) 



+ nN^ A P R x 0A lx + l a (N^ x P L X°Jx + ^{N R v) ) x A P R x° A V x 
+ d a {(C^) AX P R + (cf x A nx + u a {(C^) X A P R + {Ct ] ) x A P L _ 

+ d a (N R d) ) x A P R x 0A d x + d a (N^)i x P L x A d x + u a (N R u) ) x A P R x 0A u x + u a (N^)i x P L x° A u x 
+ h.c, (A23) 



x 



where the coefficients are 



{CI 



{l)\AX 



X 



{l)\AX 
L la 



~-92 



(rt(v)\x 

V^R laA 



(nP) ax 



(Md)\AX 
V U R la 

{C ^AX 
laA 



\f2m w cos j3 
-92(U L )A 1 (U?) X a 



(ui) A M) x a , 



(A24) 
(A25) 

(A26) 



92 

V2 



{-{Un)a 2 ~ (U N ) A l tan^} (Uf) x + - ^ a (U N ) A \U?) x 



92 m c 



y/2 \ mw cos (5 



mw cos (3 

(u* N ) A 3 (u-;) x a + 2(u* N ) A 1 t^e w (u^) 



l > a+3 



*\X 



a+3 



^= {{u N ) A 2 - {u N ) A l t^e w } (t/?) 



*\X 



-92{-{U R ) A m) X a + 



m 



=92 



\f2mw cos (3 

92{-{u L ) A \m) x + 



\f2m-w sin (5 

(Ul) A 2 (Ui) X a, 



^;(U R ) A 2 m x a+3 



\[2m w cos (3 



(u L )A 2 m x a+3 



(C 



{u)\X 
L laA 



~-92 



(u R )A 2 m x 



\^2m w sin (3 

92 



aA 



V2[ 



(U N ) A 2 + \{U N ) A X tan^} (U$ x a + ^([/,)/(^) 



(A27) 
(A28) 

(A29) 

(A30) 

(A31) 

(A32) 

(A33) 

(A34) 



x 
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{N (u ))AX 



92 

V2 



iu N ) A 2 + l -{u N ) A l t^e w \{ui) 



X 



+ 



{Us) A \Ul) x a 



raw sin (3 



+3 



4 {^^(^) A 4 (^) X Q - \{U*n)\ tan6 w (Ui) x a+3 
V2 [m^sinp 3 ^ 



(A36) 
(A37) 



APPENDIX B: DETAILS FOR e'S 

We here show the explicit form of the nLFV parameters e* T , e™ T and e^ T in the MSSMRN. 
They are calculated from 1-loop diagrams. 



1. Fore* r 



The effective coupling comprises two kinds of contribution; the one comes from the 
penguin-type diagram associated with W boson and the other is the box-type diagram: 



-y,T 



\ e u.r w 



p, T ) W -penguin i 



e ur J box- 



The penguin-part is represented as 

V c /iT/ VK-pengum /-^ t iT ' 
i 

and each contribution which is shown in Fig. |B] is calculated to be 



i 



^ ( 4vr )2 
As-3) _J K 

Pa (4vr) 2 m 2 - ml 



a(s-4) = _l I is 

Pa (4vr) 2 m 2 - ml 



{C^)i X {cT) AX D{m^M, 



x A - 



A 



(s-5) 

/3a 



1 '-t( v )\X (aA 1 )*\P A (TT \ 7=l-3/ rr h 



(An) 



X 



A 



( S -7) 

/3a 



(Bl) 
(B2) 

(B3) 
(B4) 

(B5) 

(B6) 
(B7) 



:bs) 
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X 



6 bb (O l *) BA G (m lx ,M^,M^) - {O r *) b A \f (m~ lx , M-+, M x% 



(B9) 



where the functions D, E, F and G are defined as 



D(m x ,M A ) = 



4 (1 - xax) 



{l - Ax A x + 3x 2 AX - 2x 2 AX \nx AX - 2(1 - XAx) 2 him x } 



E(m x ,m Y ,M A ) 



1 1 



x A y In xax _ x AX hi x A y 
2 x A y - x AX V 1 - a; ax 1 - a;Ay 

F(m x ,M A ,M B ) = \nx AX + 



( x\ x In x A x x 2 bx In x BX 



G(m x , Ma, M b ) = y / x Ax x BX 



xax ~ x BX \ I- xax 1 - x BX 
1 / xax hi x AX x BX In x BX 



(BIO) 
(Bll) 

(B12) 

(B13) 



^ ax - ^bx V 1 - ^ax 1 - x BX 
with xax = Mjl/m^-. The couplings for the chargino-neutralino-iy-boson interaction, L 
and R , are given asj^ 



(o L ) AB = --L(u N ) A \u R ) B 2 + (u N ) A 2 (u R ) B \ 

(O r ) a b = ^(U* N ) A 3 (U L ) B 2 + (U* N ) A 2 (U L ) B \ 



(B14) 
(B15) 



The box-part is represented as 

(<Jbox = ]T (Bie) 

i 

and each contribution, which is shown in Fig. [7J is calculated to be 

= (M^M^m^m^) (^^(ivgi^^)^^)^, (BIT) 



Bt 2) =—^jjM^,M^m^,m ffr ) (Cf ^(cf )f{N^)f (J\#>)* , (B18) 



4 v / 2G i 



" Q 4^2G 



(B19) 

l^-h (M~-,M^,m> x ,m lY ) M,-M^(C<f )^(4^)^(^ ) )^(^ ) )L, 

(B20) 



where 1^ and J4 are the functions which are given as 

J 4 (Ma, M B) mx, my) = J ^ (fc2 _ ^ )(fc2 _ M ^ {k , _ <)(fc2 _ ^ , (B21) 

/• d 4 fc k 2 
MM A ,M B ,m x ,m Y ) = J (27r)4i (Jfc2 _ M 2 )(A;2 _ M|)(F _ <)(fc2 _ m 2 y) ■ (B22) 
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X 



(5-2) 




X°A 



Vx - J * Va 



(s-5) 





lx 




FIG. 6: Dia grams which. contribut6 to (€M7-)vF-penguin 



2. For e 



Since the matter of the Earth is the neutral for U (l) e m, there is no contribution to e™ T from 
photon-penguin diagrams. The Z-penguin contribution associated with a proton and that 
with a neutron cancel each other out. The contributions which need to be taken into account 
are the Z-penguin contribution associated with an electron and the box contributions: 



\ fj/T JZ-pcnguin T I^/xt Jbox; 



,m,?t /m,u\ 

— Jbox? 



-I IT 



~ilT 



•mA ( in,d\ 

e„z = [e,j Jbox- 



~ilT 



(B23) 
(B24) 
(B25) 



The penguin contribution consists of diagrams which are drawn in Fig. |Hl 



I Z-pcnguin 



£4 



m-i) 
/it • 



(B26) 



Each diagrams is calculated to be 
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FIG. 7: Diagrams which contribute to (e*) 



box' 



A 
A 



(m-l) 


1 1 


a/3 


(4vr) 2 2 


(m-2) 


1 


a/3 


(4tt) 







\ (2sin 2 ^ - 1) (CW)^^)^^.^; 



-i + sin 2 



(B27) 



%J (^)x 7=1 " 3 (t//) 7 =i- 3 y + (sin 2 ^) {Ui)x ,=A -\u\) 1= 



Y 



=4-6 



x E(mr„,mr.,,M~ /+ 



A 



(m-3) 

'a/3 



= " T^p (2sm 2 ^ " 1) (^JJ^)? 
x 



(0" L ) A B G {m 0x ,M^, - \{0" R ) A B F (m^, il% 



-\{0' R ) A B F (m lx ,M~ +A ,M~ t 



(B28) 



(B29) 



x 



Z-boson and neutralino-neutralino-Z-boson are 



The couplings for the char gino- char gino 

(0' l ) A b = ~ {U*r)\{Ur)b 1 - \{U* R )\{U R ) B 2 + ^sin 2 ^ 
(0 ,R ) A B = - {Ul)\{U L ) B x - l -(U L ) A 2 (U L ) B 2 + SimtfOw, 
(o" L U B = - \{U N ) A \W N )\ + \{U N ) A \W N )\, 



(B30) 

(B31) 
(B32) 
(B33) 
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Xa 



(m-1) 




Xa 



(m-A) 



Z 




lx 




X A 





vx 



(m-B) 




Xa 




(m-B) 



FIG. 8: Diagrams which contributes to (e^^z-pengum- The diagrams (m-A) and (m-D) are not 
counted into e™, which are included in e* T (cf. (s-2) and (s-1) in Fig. [BJ. The diagrams (m-B) 
and (m-C) cancel each other out after the off-shell prescription (Eq. (|26|)). 



(o 



^q//L*\A 



B 



B' 



(B34) 



Here, we take into account the procedure to resolve the double counting problem which is 
explained in Sec. IIII Al 

The box contribution associated with the electron in the Earth's matter is 



-(IT 



/box 



(m,e-i) 

(IT i 



(B35) 



where 



B 



B 



(m,e-l) 

a/3 

(m,e-2) 
a/3 



^-J 4 (M,o,M,o fl ,m, x ,m zV ) (^^(JVff*)^^)^^)^, (B36) 
1 



4 v / 2G i 



(B37) 
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(B38) 



(m,e-4) T. ( AA „ A/T „ ™ . ™~ A ( Ar^h* ( At( 1 )\Y ( AT (v)*^aA ( AT (l)*^eA 

a/3 



^-J 4 (M^M^m^m^j {N%>)i B (NZyrtm } )T(K> *)y , (B39) 

(B40) 

Si7" 6) = - (M- +A ,M- p m Ix ,m, Y ) (CP)fM } )r(ClhAr(C¥y x A , (B41) 



(B42) 



w(m,e-8) 



J— J 4 (m^.Mjo .m^m^) (^)f x (JvM)J B (C7g ) *)i y (C7M*)^ J (B43) 



J a/3 



(B44) 

^7" 10) =^; J4 ( M *V M x^™->^v) (CP)f X (C^ B (Ntr Y A (nPTx A - (B45) 
The box contribution associated with the down-quark in the matter of the Earth is 



( e /i-r )box 



52Bjp+0, (B46) 

i 

where 

B^t* =^4^ (M^,M A ,m, x ,m SY ) (N^N^f (N^) Y dA (N^)t, (B47) 



8\/2Gj 
J 



(B48) 
(B49) 

47" 4) = " (M,o ,M,o ,m, x ,m dV ) (^^ fl (^)^(^>)^(JV<^)^ , (B50) 



(B51) 
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B% 4 ~ &) = - (M- p M- p m Ix ,m, Y ) (Cf >)f fl(C f ) J Y (Cf)l y (Cf)f (B52) 

The box contribution associated with the up-quark in the matter of the Earth is 



where 

^g(m,M-l) 



op "^; J4 (M,o ,M,o s ,m^,m, y ) (iV^)^(^>)f (^)^(^)^, (B54) 

(B55) 
(B56) 
(B57) 

^ =^-/* { M xP M xP m i^ m ir) (C { nYuA(C^r x A , (B58) 

(B59) 



T(m,u-6) 

5 o/3 



3. Fore£ T 

We consider the charged current interaction between neutrino beam and the nucleon in 
the detector as a detection process. It consists of the penguin contribution and the box 
contribution; 

e tr = (e^r V-penguin + (<^ T )box- (B60) 

The penguin contribution can be represented as the complex conjugate of that for e s Tfl . 
However, we must eliminate the diagrams which are already counted in the calculation of 
. The detail is shown in Sec. IIII Al 



(e^V-pengui^E^*" (B61) 
i=3 
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The box contribution is calculated to be 

(^)box = } , (B62) 

i 

where 

=^7^4 (M^M^m^m^) (cf ) B P X {Cfy BY {N^ A {N^Y x \ (B63) 
^S 2) =^T^ (M^,M,o fl) m rx ,m Jy ) (<)J B (<*)f (4" ) )^(^ ) *)x A , (B64) 

(B65) 

*C =^- p I * { M ^ M x%^l x ^r) M,-M^(iV«)^(iV^)r B (4 d >)l y (4^)f . 

(B66) 
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